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The magnetic ordering due to the long range dipole coupling in an ensemble of magnetic 
^ \ islands is investigated. If the islands are large enough and closely separated, the average 
g \ dipole energy per island can explain the magnitude of the observed ordering temperature of 
q 1 such an ensemble (U. Bovensiepen et al., J. Magn. Magn. Mater. 192, L386 (1999)). The 
energetical degeneracy with respect to a continuous in-plane rotation of the magnetic mo- 
ments in a periodic ensemble of islands is lifted in presence of an island size dispersion and an 
^ ! irregular island array. Many different (metastable) magnetic states are obtained, reminiscent 
of a spin-glass behavior. We obtain that the average magnetic binding energy per island due 
^ ' to the dipole coupling increases with increasing positional disorder. The island ensembles 
O ' exhibit non-collinear magnetic structures, resulting in non-saturated ensemble magnetiza- 
^ . tions. The calculations are performed with a classical spin model for ensembles of islands in 
! unit cells with periodic boundary conditions. The point dipole sums are augmented by an 
c3 ' island areal correction. 

s 

*^ ■ Introduction The investigation of interacting magnetic nanoparticles is a very active 
^ ■ field of current research both experimentally and theoretically |lj — ||. Measurements on 
(j growing Co/Cu(001) thin films exhibit a global ordering temperature of the order of 50 - 
£a! \ 100 K, and a magnetic hysteresis and remanence also for coverages below the percolation 
threshold Q. STM images taken in this coverage range yield an ensemble of double-layered 
£ ■ Co islands with lateral size ~ 5 nm and similar island separations, which should behave as 



a superparamagnetic ensemble. Blocking effects due to anisotropy barriers estimated by the 
Arrhenius-Neel-model |9[] could be excluded in this temperature range. In this contribution 
we will examine whether the long range magnetic dipole coupling can be responsible for a 
long range magnetic ordering at these temperatures. 

It is still under debate whether the magnetic dipole coupling can induce an equilibrium 
ordered state in an inhomogeneous two-dimensional (2D) or three-dimensional (3D) spin 
system 0-0. Certainly, a simple ferro- or antiferromagnetic state will not evolve. Instead 
a spin-glass-type |10| (random magnet) is to be expected due to the nonuniform magnetic 
interactions, which are caused by the size- and shape dispersion as well as by the positional 
disorder of the island ensemble. In such a system many different metastable states exist, 
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which are often characterized by extremely long relaxation times (magnetic viscosity). 

The influence of the dipole coupling on the blocking effects have been studied extensively 
in the Arrhenius-Neel-framework for single particles, in particular for the investigation of 
dynamical processes. The results are satisfactory for weak interactions as compared to the 
anisotropy barriers of a single island @ - . For strong interactions, however, this model is 
no longer applicable. Here we will study the strongly interacting case, considering solely the 
dipole interaction. Note that the long range indirect exchange (RKKY-) interaction might 
also be responsible for the observed magnetic ordering in ensembles of magnetic islands. 

Before describing the model let us first estimate whether the dipole coupling can account 
for the observed magnitude of the ordering temperature Tc, which must be of the order of the 
average magnetic binding energy per island. A pair of single spins has a dipole interaction 
strength of the order E dip ~ ^"i t / a o ~ 1 K, with fx at ~ 1 ji B the atomic magnetic moment, 
fiB the Bohr magneton, and a Q ~ 2.5 A the interatomic distance of 3d- transition metals. 
All iV spins of a magnetic island are aligned by the strong direct interatomic exchange cou- 
pling, and can be viewed as a single giant spin M ~ N ■ fi at (Stoner-Wohlfarth particle |11 



The diameter of a flat (2D) islands is L ~ y/N. In a densely packed island ensemble the 
island-island separation is R ~ L, resulting in an average magnetic dipole energy per island 
pair Edip ~ M 2 /R 3 ~ (fi at /al) \/N . Thus, for a typical island size of about N ~ 1000 
atoms the binding energy is Edip ~ 30 K, yielding the correct order of magnitude. In the 
following we will determine E d i P for a 2D island ensemble as functions of the size dispersion 
and positional disorder. 

Theory For the calculation of the magnetic binding energy due to the dipole coupling 
in a 2D island array we consider a unit cell with n non-overlapping, disk-shaped islands, 
and with periodic boundary conditions. Within the unit cell the islands can be placed using 
three different types of arrangements: (i) periodic square array; (ii) disturbed array, the 
island centers deviate from the sites of the periodic array using a Gaussian distribution with 
a standard deviation oy; (iii) random setup. The island size, which determines its lateral 
extension, can be either the same for all islands, of dispersed around a mean value N using 
also a Gaussian distribution with a standard deviation a^. For such an island ensemble 
we consider the following dipole-dipole Hamiltonian with classical magnetic moments Mj 
(|Mi| = N i - f i at ): 



* i,i 1 ij L 



M,,M i 4-3(r ij M l )(r ij M J 



Tij = |ry| = |rj — rj\ being the distance between the centers of gravity of islands i and j. 
The infinite range of the dipole interaction is taken into account by applying an Ewald type 
summation over all periodically arranged unit cells of the infinitely extended thin film. Since 
for a large coverage the islands are closely separated, we consider in addition to the usual 
point dipole sum the so-called areal correction (dipole-quadrupole interaction). This is of 
the order of (Ai+Aj)/rfj, where Ai cx Ni is the area of island i jl2]| , and amounts up to 50 % 
of the average dipole energy. For spherically shaped islands this correction vanishes since a 
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sphere has no quadrupole moment [fT3]] . We neglect here the lattice anisotropy. The shape 
anisotropy due to the dipole coupling within a single island vanishes for disk- or sphere- 
shaped islands. For the atomic magnetic moments and the interatomic distances we choose 
values appropriate for Fe: /i at = 2.2 /jb and a Q = 2.5 A. 

Starting from an arbitrary initial guess the total magnetic energy is minimized by varying 
the magnetic directions of each island. This magnetic relaxation is performed with the help 
of a conjugated gradient method. For simplicity, we restrict the island magnetizations to be 
always directed in-plane, Mj = Mj(cos0j, sin^j, 0), characterized by its in-plane (azimutal) 
angle fa. To account for the many different energy minima in case of a nonuniform spatial is- 
land arrangement or in presence of an island size dispersion, the magnetic energy is averaged 
over many different initial guesses for the same island arrangement. In addition, we average 
Edip over ten different realizations of the unit cell, using the same global variables (average 
sizes, standard deviations, etc.) characterizing the island ensemble. The energy reference 
is given by a random set of magnetic angles {fa}, refering to a completely disordered system. 



Results Before presenting the average magnetic dipole energy in particular for a disordered 
island array, let us recall shortly the resulting structure for a uniform system, i.e. the islands 
having all the same size and placed on a square mesh. Two different magnetic arrangements 
are possible: (i) a metastable parallel magnetization of the islands (ferromagnetic solution), 
and (ii) a columnar arrangement consisting of magnetized rows (or columns) of islands with 
alternating magnetic orientations. The latter is the ground state of a square array of spins 
or magnetic islands, having a vanishing remanent magnetization. Note that both solutions 
exhibit a continuous degeneracy with respect to an in-plane rotation. This property is 
one of the prerequisites for a vanishing global magnetic order of a 2D magnetic system at 
finite temperatures (Mermin- Wagner-theorem [fLlp . By introducing now a size- or positional 
disorder into the island ensemble, this degeneracy is immediately lifted. The system exhibits 
a (possibly very large) number of discrete states, which are separated by energy barriers 
and have approximately the same energy if no external magnetic fields or aligned anisotropy 
easy axes are present. The magnetic arrangement of the island ensemble is strongly non- 
collinear, the net magnetization is small or vanishes completely (demagnetizing effect due 
to flux closure), cf. Fig.l. We emphasize that these non-collinear arrangements do not 
correspond to a disordered system, since the local magnetizations are strongly correlated by 
the dipole interaction. The parallel and the columnar magnetic arrangements of the island 
ensemble are always unfavorable, and even do no longer refer to a stable state. Furthermore, 
due to the lifting of the energy degeneracy the above mentioned prerequisite for the Mermin- 
Wagner-theorem is no longer fulfilled. Therefore, a dipole induced magnetic order might exist 
at finite temperatures also for an inhomogeneous 2D system ('order- by-disorder-effect' [p"5|j). 

In Fig. 2 we show the average magnetic dipole energy E^p per island as a function of the 
island positional disorder, characterized by the standard deviation o r . The (average) island 
size is chosen to be N = 100, the average island-island distance is 1.5 times the (average) 
island diameter, resulting in a coverage of about 35 %. E<n v decreases with increasing 
positional disorder, reaching its minimal value (= maximal magnetic binding energy per 
island) for a random setup of the island ensemble. The increase of binding energy with 
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respect to the uniform system amounts to roughly 50 %. As can be also seen from Fig.2, an 
island-size dispersion has no large effect on E^ip- 

The reason for the different effects of the two types of disorder on the magnetic binding 
energy is that E d i P is a bilinear function of the island sizes iVj, the size dispersion effect 
averages out for symmetrically distributed island sizes around the mean value N. On the 
other hand, the dipole energy has a nonlinear dependence on the island-island distance: 
E dip oc r~j 3 . Thus, with increasing positional disorder the decrease of \E dip \ with an enlarged 
distance between some island pairs is more than counterbalanced by a corresponding 
increase of \Edi P \ for smaller distances between other island pairs. This leads to an increase 
of the average magnetic binding energy of the island ensemble. 



Conclusion We have calculated the magnetic binding energy due to the long range dipole 
coupling in a 2D ensemble of magnetic islands. The binding energy increases by ca. 50 % for 
a random arrangement of islands with respect to a periodic array. An island size dispersion 
has less effect on the binding energy. Due to the nonuniform system a large number of 
different metastable states emerges with strongly non-collinear magnetic arrangements, as 
reminiscent of a spin-glass ||10|| . We expect that if the islands are large enough and closely 
separated, the dipole energy can explain the observed magnetic ordering of a growing thin 
film below its percolation threshold As a first approximation, the resulting ordering 
temperature can be determined from a mean-field approximation, and corresponds to the 
spin-glass temperature. For a comparison with experiments also magnetic anisotropies need 
to be taken into account, as well as overlapping magnetic islands [If]]. Furthermore, one can 
calculate the ensemble magnetization and susceptibility as a function of an applied magnetic 
field. 
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Figure 1: Visualization of a typical (meta-) stable magnetic arrangement in the unit cell 
as obtained from the calculations for the average dipole energy. The average island size is 
N = 100 spins, yielding an average island radius of about 10 a , a the interatomic distance. 
The size of the unit cell is chosen to be (150 a ) x (150 a ), containing 100 islands, and 
yielding a coverage ~ 35 %. The spatial island arrangement corresponds to a disturbed 
array, the island centers deviate from the sites of a square array with a standard deviation 
a r = 10 %. The island size exhibits a dispersion of about o^r = 20 %. The arrows depict 
the in-plane directions of the island magnetizations. Note the tendency to the formation of 
island chains and flux closure structures. 
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Figure 2: Average dipole energy E dip (a r ) per island of a 2D ensemble of magnetic islands as 
a function of the positional disorder, which is measured by its standard deviation a r from 
a square array, s. Fig.l. The full line refers to a system with the same size for all islands, 
the dashed line to a system with 20 % size dispersion. The line denoted by 'random' refers 
to Edip of a random island arrangement in the unit cell. For the atomic magnetic moment 
and the interatomic distance we have choosen values appropriate for Fe: /i at = 2.2 /ib and 
a Q = 2.5 A. 
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